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j—{ , Abstract 
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^SJ ' Within a perturbative cosmological regime of loop quantum gravity corrections to 

. effective constraints are computed. This takes into account aU inhomogeneous degrees 

I of freedom relevant for scalar metric modes around flat space and results in explicit 

' expressions for modified coefficients and of higher order terms. It also ihustrates the 

, role of different scales determining the relative magnitude of corrections. Our results 

[ demonstrate that loop quantum gravity has the correct classical limit, at least in its 
sector of cosmological perturbations around flat space, in the sense of perturbative 

^ ■ effective theory. 

bp: 

1 Introduction 

rS ■ 

' Interacting quantum theories in low energy or semiclassical regimes can be described by ef- 
fective equations which amend the classical ones by correction terms. Compared to the full 
quantum description an analysis of such effective equations is much easier once they have 
been derived. In addition to the simpler mathematical structure, difficult conceptual and 
interpretational problems of wave functions can be evaded, still allowing one to compute 
potentially observable effects. Technical and conceptual problems are even more severe in 
quantum gravity, in particular in background independent formulations. Yet, especially in 
this case observational guidance would be of invaluable help. Since the high energy regimes 
of cosmology are commonly considered as the best access to such guidance, an effective 
description for fields relevant for early universe cosmology is needed. In this paper we 
use an effective framework of perturbative loop quantum gravity around a spatially flat 
isotropic background space-time to derive correction terms to the classical constraints. 
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Our analysis will be done for scalar metric modes in longitudinal gauge as this can be 
used to simplify the perturbative basic variables. They can then be chosen to be of diagonal 
form, although now fully inhomogeneous, which is the main reason for simplifications as 
they have been used extensively in symmetric models [H [2], [3l H] . The main constructions 
of these models can thus be extended, in a similarly explicit form, to inhomogeneous 
situations without assuming any symmetry. This allows us to compute explicit corrections 
to effective constraints which, in combination with the Hamiltonian analysis of cosmological 
perturbation theory in [5], leads to corrected perturbation equations and new effects f6j. 

A physical regime is selected by introducing a background geometry in the background 
independent quantization through a specific class of states [7j . This keeps the characteristic 
background independent features of the quantum theory, such as its spatial discreteness, 
intact while bringing the theory to a form suitable for perturbative expansions and appli- 
cations. In the perturbative regime, we will make use of special structures provided by 
the geometrical background which can usually be chosen to allow symmetries, e.g. isotropy 
for cosmological perturbations around a Friedmann-Robertson- Walker model. In particu- 
lar, we use this to introduce regular lattice states with a spacing (in geometrical variables 
rather than embedding coordinates) whose size determines at which scales quantum effects 
become important. The geometrical spacing thus specifies on which scales physical fields 
are being probed by a given class of states. On such a lattice, explicit calculations can be 
done. 

We demonstrate this by providing higher curvature terms as well as corrections to 
inverse powers of metric components. Several issues that arose in isotropic models will be 
clarified. Finally, we discuss general aspects of effective equations and the semiclassical 
limit of loop quantum gravity. The article thus consists of two parts, an explicit scheme to 
derive correction terms presented in Sec. [3] and IH and a more general discussion of effective 
equations and the classical hmit in Sec. O 

2 Basic variables and operators 

The basic variables of interest for a canonical formulation of gravity [S] are the spatial 
metric qab occurring in the space-time metric 



(or equivalent quantities such as a triad e\ with e^e^ = g^b oi its inverse ef), extrinsic 
curvature Kah = {'^N)~^{qab — DaN^ — D^Na) (or related objects such as the Ashtekar 
connection) and matter fields with their momenta. The components (lapse function) 
and N°- (shift vector) of the space-time metric are not dynamical, and thus do not have 
momenta, but are important for selecting the space-time slicing or the gauge. 

Because of their role in loop quantum gravity, we will use Ashtekar variables [9l [TO] 
which are a densitized triad E°- = \ det el\e^ and the connection = — •jKl^ with the 
spin connection 



ds^ = -N^t^ + Qabidx" + N''dt){dx'' + N'^dt) . 



(1) 




(2) 
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compatible with the triad, K^^ = Kabe\ and the positive Barbero-Immirzi parameter 7 [TOl 
[TT] . We use them here in perturbative form on a flat isotropic metric background, focusing 
on scalar modes. This means, as explained in more detail in [5], that the unperturbed 
metric as well as its perturbations can be assumed to be diagonal, = p^'^\x)61, which 
simplifies calculations. For scalar modes, all diagonal components of the metric qab = 
a^(l — 2ip{xY)6ab are in fact equal, but we will see that this restriction is not general 
enough for formulating a loop quantization. Moreover, we can choose a vanishing shift 
vector A^"^ = 0, implying that extrinsic curvature K^^ = k{i){x)5l^ is diagonal, too. (The 
Ashtekar connection, on the other hand, will not be diagonal because it has non-diagonal 
contributions from the spin connection. It is of the form = k(^i){x)5\ + ipi{x)e'^J where 
the non-diagonal part ipj arising from the spin connection computed in Sec. 13.1.41 can be 
dealt with perturbatively.) Our calculations will thus be done in a given gauge, and would 
be more complicated in others. Nevertheless, we are including the general perturbations 
of metric and matter variables relevant for scalar modes without too strong restrictions as 
they could arise in other gauges. 

2.1 Gauge choices and their implications for a quantization 

In general, the space-time gauge is determined by prescribing the behavior of lapse function 
N and shift vector iV" occurring in a metric ([H). The lapse function, as we will see, can 
be chosen arbitrarily in our calculations, but the shift vector is restricted for a diagonal 
perturbation to be realized. We are thus using a particular class of gauges in setting up our 
calculations, although we do not explicitly make use of the form of gauge transformations. 
This is important because the canonical constraints, most importantly the Hamiltonian 
constraint H in addition to the diffeomorphism constraint Da, and thus also the gauge 
transformations 6^f = {f,C,^H + ^"-Da} they generate will be corrected by quantum ef- 
fects. Classical properties of the gauge transformations should thus not be used before one 
computes quantum corrections. It is then a priori unclear which particular gauge choices, 
other than fixing lapse and shift directly, are allowed. Some gauges implicitly refer to 
gauge transformation equations to relate metric to matter perturbations, or to select the 
space-time slicing such as for the fiat gauge. In this example, one would make use of gauge 
transformations to set the spatial metric perturbation equal to zero which allows one to 
focus calculations on the simpler matter part. In this process, one solves gauge transfor- 
mation equations of the metric perturbation, depending on lapse and shift, such that the 
transformed perturbation vanishes. This determines a gauge to be chosen, but makes use 
of explicit gauge transformation equations which are not guaranteed to remain unchanged 
with quantum corrections. Our choice of vanishing shift, on the other hand, is harmless 
because it does not refer to explicit gauge transformation equations. We are thus working 
at a more general level keeping metric and matter perturbations independent. A combined 
gauge invariant combination of the two perturbations can be determined once the quantum 
corrected gauge transformations have been computed. 

When constraints are modified, manifest covariance of the resulting equation becomes 
an issue as it is discussed in more detail in [5] . Such quantum corrections are derived from 
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effective constraints of gravity which are defined as expectation values of quantum gravity 
operators in general states [12j. We motivate the procedure here briefly, and will provide 
some further information in Sec. [SI for details we refer to [T2 l [T3 l [T4] . If constraint operators 
satisfy the classical constraint algebra, covariance would be manifest for the expectation 
values. But there is an additional step involved in deriving effective equations: the expec- 
tation values depend on infinitely many quantum variables such as the spreads of states 
which do not have classical analogs. Effective equations are obtained by truncating this to 
finitely many fields (similarly to the derivative expansion in low energy effective actions), 
resulting in equations of motion of the classical form corrected by quantum terms. Indeed, 
any quantum theory is based on states which are not just determined by expectation values 
of the basic variables such as and E'f in loop quantum gravity. Expectation values of 
the basic variables would correspond to the classical values in constraint expressions, but 
there are infinitely many further parameters such as the spread and deformations of the 
state from a Gaussian. These additional, quantum variables can suitably be parameterized 
in the form 

= ((^-(^)r"(p.-(p.)r)weyi (s) 

for any degree of freedom {q,Pq) present in the classical theory. Here, l<nGN, 0<a<n, 
and the subscript "Weyl" denotes totally symmetric ordering. Every classical degree of 
freedom thus does not only give rise to expectation values (g) and (pg) but to infinitely 
many additional quantum variables. All of these variables are dynamical, and are in general 
coupled to each other. 

Moreover, expectation values of most operators, including Hamiltonians, in general 
states depend on all these infinitely many variables. This is to be reduced to a finite set 
for an effective description which introduces quantum correction terms into the classical 
equations. In particular, spread and deformation parameters are usually assumed to be 
subdominant compared to expectation values. Without explicitly constructing semiclassi- 
cal states satisfying such conditions, one can make semiclassicality assumptions for those 
parameters to be negligible. This is what we will do in this paper as a shortcut to deriving 
effective expressions from a full quantum theory. Since special quantization steps are in- 
volved in the construction of operators which reformulate classical expressions, corrections 
in effective equations will result which are not sensitive to the precise form of semiclassical 
states. 

2.2 Lattice states and basic operators 

We are thus able to implement all degrees of freedom needed for inhomogeneities in a 
way which is accessible to explicit calculations. While the general kinematical arena of 
loop quantum gravity is based on discrete spatial structures built on arbitrary graphs with 
possibly high-valent vertices, we will use regular lattices with 6-valent vertices. Regularity 
of the lattice is implemented by making use of symmetries of the background we are 
perturbing around: The three independent generators of translational symmetry define 
lattice directions. In explicit constructions of lattice states, a scale Eq will appear which 
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is the coordinate length of lattice links measured in a given, fixed embeddingQ But this 
parameter is independent of the quantum variables assigned to each link we will be using, 
which means that the quantum theory will be defined on "freely fioating" lattices as in the 
full theory, respecting diffeomorphism invariance. The scale will only become important 
in the continuum limit, when making contact between the quantum variables and classical 
continuous fields. This obviously breaks manifest diffeomorphism covariance, just as the 
classical perturbation theory in basic fields rather than gauge-invariant combinations, since 
the classical perturbations are written with respect to a background space-time. 

Compared to the full theory, we are restricting states by assuming regularity and thus 
allowing, e.g., only unknotted links and vertices of valence at most six. This turns out 
to be sufficient to include all relevant perturbative degrees of freedom. While the general 
graphs of loop quantum gravity allow more freedom, its meaning is not known and appears 
redundant in our application. 

2.2.1 Holonomies and fluxes 

The canonical fields are given by (A^, E^) which are to be turned into operators on a suitable 
Hilbert space. To set this up, we need to choose a functional representation of state, which 
is conveniently done in the connection representation where states are functionals of A\. 
According to loop quantum gravity, lattice graphs then label states and determine their 
expressions as functions of connections: a state associated with a given graph depends on 
the connection only through holonomies 

=Pexp(/dtA^eV,) 

e 

along its edges. Here Tj = — are the S'?7(2)-generators in terms of Pauli matrices aj 
and V denotes path ordering. That those are the basic objects represented on a Hilbert 
space together with fiuxes 

Fs{E) = [ d'yE^r'na (4) 

for surfaces 5* with co-normal is the basic assumption of loop quantum gravity [15j. Our 
corrections to cosmological perturbation equations will be implications of this fact and thus 
test the theory directly. Using the perturbative form of A^, we can split off perturbatively 
the non-diagonal part (composed of spin connection components) in an expansion and 
exploit the diagonality of the remaining part to obtain h^j = exp(7r/ ^dtkj{ejjj(t)). 

Similarly, fiuxes will be of the form Fyj = Js ^ ^yp\y)- A lattice link starting at a vertex 
V in direction / in a fixed orientation is denoted e^, /, and a lattice plaquette transversal to 
this edge and centered at its midpoint as S^j. (Here and in the following set-up we closely 
follow [7] to which we refer the reader for more details; note, however, that 7 has been 
absorbed there in fc/.) 

^The coordinate length need not be the same for all links, but can be chosen this way without loss of 
generality. 



5 



Matrix elements of the variables h^j together with Fyj form the basic objects of loop 
quantum gravity in this setting. They are thus elementary degrees of freedom, comparable 
to atoms in condensed matter. Classical fields will, as we display in detail later, emerge 
from these objects in suitable regimes and limits only. Even in such regimes where one can 
recover the usual metric perturbations there will in general be correction terms examples 
of which we aim to compute below. Correction functions will then also depend on the 
basic objects h^j and Fyj directly, which can be expressed through the classical metric 
perturbations in a secondary stepo 

To recover the correct semiclassical behavior one has to make sure that effective equa- 
tions of motion can indeed be written in a form close to the classical ones. Since classical 
Hamiltonians are local functionals of extrinsic curvature and densitized triad components, 
it must then be possible to approximate the non-local, integrated objects h^j and Fyj by 
local values of kj and evaluated in single points This is indeed possible if we assume that 
kj is approximately constant along any edge, whose coordinate lengths are io = f dt. 
We can then write 

Kj = exp( / dt-fky) ^ cos{io-fki{v + J/2) /2) + 2tj sm{£o-fki{v + 1/2) /2) (5) 

where v + I /2 denotes, in a slight abuse of notation, the midpoint of the edge which we 
use as the most symmetric relation between holonomies and continuous fields, and 

F,j= j p\y)d^y^tlp\v + I/2) (6) 

(note that the surface S^j is defined to be centered at the midpoint of the edge e^j). 
This requires the lattice to be fine enough, which will be true in regimes where fields 
are not strongly varying. For more general regimes this assumption has to be dropped 
and non-local objects appear even in effective approximations since a function ki will 
be underdetermined in terms of the h^j. Since the recovered classical fields must be 
continuous, this means that they can arise only if quantizations of h^j and F^j, respectively, 
for nearby lattice links do not have too much differing expectation values in a semiclassical 
state. If this is not satisfied, continuous classical fields can only be recovered after a process 
of coarse graining as we will briefly discuss in Sec. 15.2. 1[ 

In addition to the assumption of slowly varying fields on the lattice scale, we have also 
made use of the diagonality of extrinsic curvature which allows us to evaluate the holonomy 
in a simple way without taking care of the factor ordering of su(2)-values along the path. 
We can thus re-formulate the theory in terms of U(l)-holonomies 

r/,,, = exp(^ / dt-iki/2) ^ exp(^47^/(^ + //2)/2) (7) 

Ed, I 

^This already indicates that difficulties which were sometimes perceived in isotropic models, where 
corrections seemed to depend on the scale factor whose total scale is undetermined, do not occur in this 
inhomogeneous setup. 
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along all lattice links e^j. On the lattice, a basis of all possible states is then given by 
specifying an integer label for each edge starting at vertex v in direction I and defining 

{k{x)\ . . . ,fiyj, . . .) ■.= Y\_Gxp{ifiyj J dt-fki/2) (8) 

v,i 

as the functional form of the state | . . . , Hv,i, ■ ■ ■) in the ^-representation. The form of the 
states is a consequence of the representation of holonomies. States are functions of U(l)- 
holonomies, and any such function can be expanded in terms of irreducible representations 
which for U(l) are just integer powers. This would be more complicated if we allowed 
all possible, also non-diagonal, curvature components as one is doing in the full theory. 
In such a case, one would not be able to reduce the original SU(2)-holonomies to simple 
phase factors and more complicated multiplication rules would have to be considered. In 
particular, one would have to make sure that matrix elements of holonomies are multiplied 
with each other in such a way that functions invariant under SU(2)-gauge rotations result 
^6j. This requires additional vertex labels which we do not need in the perturbative 
situation. 

For the same reason we have simple multiplication operators given by holonomies as- 
sociated with lattice links, 

Vvjl f^v',J, . . .) = I . . . , + 1, . . .) . (9) 

There are also derivative operators with respect to kj, quantizing the conjugate triad 
components. Just as holonomies are obtained by integrating the connection or extrinsic 
curvature, densitized triad components are integrated on surfaces, then called fluxes (jlj), 
before they can be quantized. For a surface S of lattice plaquette size intersecting a single 
edge e^j outside a vertex, we have 

Fv,i\---, IJ'v',j, ■■■) = 47r74/Ui,,/| • • • , fJ'v^j, ■ ■ ■) (10) 

or 

^v,i\ ■ ■ ■,(J'v',j, ■■■) = 27r74(/it,-/,/ + IJ'v,i) \ ■ ■ ■,fJ'v',j, ■■■) ■ (11) 

if the intersection happens to be at the vertex. The Planck length ip = \jGfh arises 
through a combination of G from the basic Poisson brackets and h from a quantization 
of momenta as derivative operators. Here, in a similar notation as above, v — I denotes 
the vertex preceding v along direction / in the given orientation. We will later call such 
labels simply ^y-ij = fiv-i as illustrated in Fig. [H These operators quantize integrated 
triad components ([6]). This shows that all basic degrees of freedom relevant for us can be 
implemented without having to use the more involved SU(2)-formulation. 

Note, that even for scalar perturbations which classically have triads proportional to 
the identity, distinct p^(f )-components have to be treated differently at the quantum level. 
One cannot assume all edge labels around any given vertex to be identical while still 
allowing inhomogeneity. Moreover, operators require local edge holonomies which change 
one edge label fi^j independently of the others. Similarly, corresponding operators Fyj 



7 



-1 

V —I ^v-I,I—^,-I ^v,I V^l 

Mv-/,/=Mv,-/ Mr,/ 

Figure 1: Edges adjacent to a vertex f in a given direction /. For the edge oriented 
oppositely to the chosen one for direction J, the labels "f — I, I" and "f , — /" can be 
chosen interchangeably, defining in this way negative values for the label /. 

and Fy^j (/ 7^ J) act on different links coming out of a vertex v and have thus independent 
eigenvalues in general. To pick a regime of scalar modes, one will choose a state whose 
edge fluxes are peaked close to the same triad value in all directions and whose holonomies 
are peaked close to the same exponentiated extrinsic curvature values, thus giving effective 
equations for a single scalar mode function. But this to equal values in different directions 
cannot be done at the level of operators. 

These basic operators h^j and F^j will appear in more complicated ones and in partic- 
ular in the constraints. As we can see, they depend not directly on the classical flelds p\x) 
and kj{x) but, in local approximations, on quantities p\x) := C-lP^x) and := ioki{x) 

rescaled by factors of the lattice link size io. This re-scaling occurring automatically in our 
basic variables has two advantages: It makes the basic variables independent of coordinates 
and provides them unambiguously with dimensions of length squared for p while k becomes 
dimensionless. (Otherwise, one could choose to put dimensions in coordinates or in metric 
components which would make arguments for the expected relevance of quantum correc- 
tions more complicated.) This also happens in homogeneous models [T7], but in that case, 
especially in spatially flat models, there was sometimes confusion about the meaning of 
even the re-scaled variables. This is because the scale factor, for instance, as the isotropic 
analog of could be multiplied by an arbitrary constant and thus the total scale would 
have no meaning even when multiplied by the analog of i^. Thus, correction functions 
depending on this quantity in an isotropic model require an additional assumption on how 
the total scale is flxed. 

This is not necessary in inhomogeneous situations. Here, the quantities will appear in 
quantum corrections and their values determine unambiguously when corrections become 
important. The corresponding fluxes are the relevant quantum excitations, and when 
they are close to the Planck scale quantum corrections will unambiguously become large. 
On the other hand, if the p^ become too large, approaching the Hubble length squared 
or a typical wave length squared, discreteness effects become noticeable even in usual 
physics. As we will see in more detail in Sec. 15.2.21 this allows one to estimate orders of 
magnitudes of corrections to be expected even without detailed calculations [B]. Although 
the size of the p^ is coordinate independent, unlike the value of the scale factor, say, its 
relation to the classical fleld depends on io and thus on the lattice size. It may thus 
appear that p^ is coordinate dependent, but this is clearly not the case because it derives 
directly from a coordinate independent flux. The lattice values are deflned independently 
of coordinates, just by attaching labels /i^, / to lattice links. Once they have been specifled 
and the lattice has been embedded in a spatial manifold, their relation to classical metric 
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fields can be determined. It is, of course, tlie classical fields such as metric components 
which depend on the coordinate choice when they are tensorial. The relation between 
and the classical metric depends on the lattice spacing measured in coordinates because 
the representation of the classical metric itself depends on which coordinates have been 
chosen. Thus, our basic quantities are coordinate independent and coordinates enter only 
when classical descriptions are recovered in a semiclassical limit. 

2.2.2 Volume 

An important ingredient to construct constraints is the volume operator. Using the classical 



expression V = /d^xv^ 

p2p3| introduce the volume operator V = n/=i y 1-^^, 
which, using ( ITTl) . has eigenvalues 



V 1=1 

While densitized triad components are directly implemented through basic fluxes, the 
process of quantizing triad or co-triad components is more indirect. While they are uniquely 
determined from the densitized triad classically, one needs to take inverse components. 
With flux operators having discrete spectra containing zero, this is not possible in a direct 
manner at the quantum level. Nevertheless [18], one can construct operators for co-triad 
components based on the classical identity 

I A^, J v^dit^ld^xj = 2n-fGe'^''eabc-j=^= = ^rc-fGel . (13) 

On the left hand side, no inverse appears and we just need to express connection compo- 
nents in terms of holonomies, use the volume operator and replace the Poisson bracket by 
a commutator. Resulting operators are then of the form he[h~^,V] for SU(2)-holonomies 
along suitable edges e, e.g. 

tT{T'Kj[h-}, K]) ~ -Umo{a[X} (14) 

for hyj as in (l5l). This shows that factors of the link size io are needed in reformulating 
Poisson brackets through commutators with holonomies, which, as will become clear below, 
are provided by the discretized integration measure in spatial integrations such as they 
occur in the Hamiltonian constraint. 



3 Hamiltonian constraint 

Holonomies, the volume operator and commutators between them are finally used to define 
Hamiltonian constraint operators. We will briefly describe the general procedure and then 
derive resulting correction terms in effective equations for both gravitational and matter 
contributions to the constraint. 
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3.1 Gravitational part 

The gravitational contribution to the Hamiltonian constraint is given by 

^f^] = Y^^d='a;iV|detEr^/^(e.,,F^,E;E^ (15) 
- 2(1 + - - n)EfE'}) 

in terms of Ashtekar variables with the curvature F^j, = 2d[aA\^ + e^^^A^^A^. The second 
term, quadratic in extrinsic curvature components Kl^ = 7~^(r^ — A^), is in general more 
complicated to deal with due to the appearance of spin connection components as func- 
tional of through ([2]). One usually starts with quantizing the first term and then uses 
the identity [18] 




which allows one to express the second contribution in terms of the first. In the first term, 
then, the densitized triad components including the inverse determinant can be quantized 
using (ITS]) , and the curvature components can be obtained through holonomies around 
appropriately chosen small loops [19]. On our regular lattices, natural loops based at a given 
vertex are provided by the adjacent lattice plaquettes. After replacing the Poisson brackets 
by commutators, the resulting first part of the Hamiltonian operator, H'^'^^ = H^^^ has 
non-zero action only in vertices of a lattice state, each contribution being of the form 

= 16^8^^^ ^ "^"^"^^ ^^^^ 

' IJKcrre{±l} 

{A)-^K^„jj (A) ^h^^aKK{.A)[h^^„j^K 

summed over all non-planar triples of edges in all possible orientations. (There are 48 terms 
in the sum, but we need to divide only by 8 since a factor of six arises in the contraction 
of basic fields occurring in the constraint.) The combination 

(A) K,ajj{A) 

gives a single plaquette holonomy with tangent vectors e^^o-j/ and e^^a-jj as illustrated in 

Fig. m 

When expanded in io assuming sufficiently small edges, the leading term is of the order 
^0 which automatically results in a Riemann sum representation of the first term in (fT5il . 
This justifies H^^^ as a quantization of the classical expression. As seen from the argument, 
one needs to assume that the lattice is sufficiently fine for classical values of the fields 
Al^. Thus, there are states corresponding to coarser lattices on which stronger quantum 
corrections can result. As usually, semiclassical behavior is not realized on all states but 
only for a select class. For any low-curvature classical configuration, one can make sure 
that a chosen lattice leads only to small quantum corrections such that sufficiently many 
semiclassical states exist. 
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Figure 2: Elementary lattice plaquette with holonomies around a closed loop. 



3.1.1 Quantization 

The required calculations for SU(2) holonomies and their products usually do not allow 
explicit diagonalizations of operators. But some physical regimes allow one to decouple the 
matrix components at least approximately. This is realized for several symmetric models 
and also for perturbations at least of some metric modes around them. In particular, 
after splitting off the non-diagonal part of the connection in the perturbative expansion 
considered here, we can take the trace explicitly and reduce the expression to U(l). Since 
the diagonal part of the Ashtekar connection for our perturbations is contributed entirely by 
extrinsic curvature, we are effectively using "holonomies" computed for extrinsic curvature 
rather than the Ashtekar connection. Although extrinsic curvature is a tensor rather than 
a connection, it is meaningful to use it in expressions resembling holonomies, denoted here 
simply as h^j, on a given metric background. This has the additional advantage of easily 
combining the remaining quadratic terms in K^^ with the first term of the constraint f[T^ 
without using squares of multiple commutators from quantizing (fT6l) . Writing 

= 2d^,Tl^ + 2jd[aKl^ + 7e.,fc {TiK^ + T^Ki) + e,,, (F^r^ + j'KiK^) (18) 

we obtain a term 2'yd[aK^ + ■j^KlK^ resembling "curvature" F^^^^K) as computed from 
extrinsic curvature alone, a curvature term of the spin connection as well as cross-terms 
^ijk ij^a^b ~^ ^b^a) ■ o^^^ coutcxt duc to the diagouality conditions the cross-terms 
disappear [5] and we only have the i^'-curvature term and spin connection terms to quantize. 
The first term can then be combined with the term quadratic in /T* in f|T5|) . removing the 
need to use double commutators. We denote this contribution to the constraint as 

Hk[N] := ^ / <^'xN |det E\-"^ {e^.kldaKl - KiK^) E^e'^ (19) 
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(since also daK^ drops out as used later, the constraint is 7-independent) and the remaining 
term as 

Ht[N] := ^ / d'xN |det E\-'/' {e^.^daTl + ^^t) E^e'^ . (20) 

Then, H[N] = Hk[N] + H-p[N] is the constraint for scalar modes in longitudinal gauge. 
Both terms can rather easily be dealt with, using holonomies around a loop for the first 
term (this subsection) and direct quantizations of for the second (Sec. 13.1.41) . The 
split-off spin connection components are thus quantized separately, which is possible in the 
perturbative treatment on a background, and then added on to the operator. 

Note also that as a further simplification the derivative term of extrinsic curvature 
disappears from the constraint for diagonal variables as assumed here. This will automat- 
ically happen also from holonomies around loops. We emphasize that the quantization 
procedure followed is special to the given context of scalar perturbations on a flat isotropic 
background. Nevertheless, it mimics essential steps of the full constructions as discussed 
in more detail in Sec. 14.2. 1[ Its main advantage is that it allows explicit derivations of all 
necessary terms and thus explicit effective equations to be confronted with observations. 
Moreover, it is far from clear that the constructions currently done in the full setting will 
remain unchanged with further developments. We thus evaluate the key features of the 
scheme without paying too close attention to current details. 

Following the general procedure, we thus obtain vertex contributions 

M 1 2i N{v) ^ ^ ^ ^ ^ JJK 

' IJKaie{±l} 



As before, h^j denotes a i^-holonomy along the edge oriented in the positive /-direction 
and starting at vertex t>, but we also include the opposite direction h^-i in the sum to 
ensure rotational invariance. Note that following our convention, such holonomies are 
identified with j. In some of the holonomies, f + / is again the vertex adjacent to v 
in the positive /-direction. The {/J/Tj-summation is taken over all possible orientations 
of the / J-loop and a transversal /T-direction. Also, for notational brevity, we introduce 

Cyj ■.= ]^ii{Kj) , Syj := -ii{T(^i)Kj) (22) 

such that (j5]) becomes h^j = c^j + 2r/s^,/. In a continuum approximation, we have 

c^,7 = cos(7A;/(t; + //2)/2) , s^j = sm{-fki{v + I /2) /2) (23) 

where kj(v) = iQki{v). After substituting this expression into fl2T]) and making use of the 
identit}o (for some fixed /, J, K and numbers Xj and yi) 

e/jii- tr [(xil + 2yiTi){x2l + 2t/2Tj)(x3l + 2y3TK)] = e/j/^ tr(xiX2X3l) + SeuK tT{yiy2y3TiTjTK) 

= 2{xiX2X3 - yiy2y3(^IJK)(^IJK, 



^Here the fundamental representation of tj has been used: tr(I) = 2, tr(T/T/) = —-^Sij 
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any one term of the sum in (12T|) becomes 

^— ^ ti{K jK+i^jh-^j^jh-^jK^KK]^, K]) (24) 

where 

Bv,K '■= ^-j^l^Qfi {^^'",KVvCy^K — C-u^KVvSy^K^ . (25) 



In the first hne of fl24l) . the expression inside the curly braces is symmetric in the in- 
dices I and J, hence vanishes when contracted with tuK- Therefore only the second line 
contributes, and the extrinsic curvature part of the gravitational constraint is 



-N{v) 

1 

} (26) 



fr _ -iVW \^ 5\( - \ 

'^'^ IJKai(i{±l} 



_ -N{v) ^ ^ _^ IR 1 

where in the last line trigonometric identities have been used to express products of sines 
and cosines through 

<a,/,<x,,j := sin {^-^{K,i{y + aill2) ± K^i{v + ojJ + oilj^) 

As in this expression, all functions ki are, as before, evaluated at the midpoint of the edge 
e„^/. We see that in the homogeneous case the first term in the sum vanishes and the 
leading contribution is 

4 sin(7fc,/2) cos(7fc,/2) sin(7fcj/2) cos(7fcj/2)5„,^, (27) 

in agreement with [2]. 



3.1.2 Higher curvature corrections 

There are two types of corrections visible from this expression: Using commutators to 
quantize inverse densitized triad components implies eigenvalues of B^^i which differ from 
the classical expectation at small labels /it,,/. Moreover, using holonomies contributes higher 
order terms in extrinsic curvature together with higher spatial derivatives when sines and 
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cosines are expanded in small curvature regimes. We will now compute the next-leading 
terms of higher powers and spatial derivatives of kj{v) before dealing with inverse power 
corrections in the following subsection. 

First, recall the usual expectation that quantum gravity gives rise to low energy effective 
actions with higher curvature terms such as / d^x^/l detglipR"^ or J d^x-y/l det g\£pRf^upcrR^'^'^'^ 
added to the Einstein-Hilbert action J d^x^/l det g\R. Irrespective of details of numerical 
coefficients, there are two key aspects: The Planck length ip = \J GH must be involved for 
dimensional reasons in the absence of any other length scale, and higher spatial as well 
as time derivatives arise with higher powers of R^upa- In canonical variables, one expects 
higher powers and higher spatial derivatives of extrinsic curvature and the triad, together 
with components of the inverse metric necessary to define scalar quantities from higher 
curvature powers (which forces one to raise indices on the Riemann tensor, for instance). 
Higher time derivatives, on the other hand, are more difficult to see in a canonical treat- 
ment and correspond to the presence of independent quantum variables without classical 
analog [13]. 

Any quantization such as that followed here starts from the purely classical action where 
h and thus ip vanishes. In effective equations of the resulting quantum theory, quantum 
corrections depending on h will nevertheless emerge. As a first step in deriving such effective 
equations, we have non-local holonomy terms in a Hamiltonian operator which through its 
expectation values in semiclassical states will give rise to similar contributions of the same 
functional form of ki{v). At first sight, however, the expressions above do not agree with 
expectations from higher curvature actions: One can easily see that in fl26l) there are higher 
powers of extrinsic curvature by expanding the trigonometric functions, and higher spatial 
derivatives of extrinsic curvature by Taylor expanding the discrete displacement involved, 
e.g., in ki{v + I /2). Moreover, higher spatial derivatives of the triad arise from similar 
non-local terms in the spin connection contribution Hp discussed later. But there are no 
factors of the Planck length in such higher powers (all factors of G and h are written out 
explicitly and not "set equal to one"). In fact, by definition ki{v) is dimensionless since 
it is obtained by multiplying the curvature component ki{v) with io in which all possible 
dimensions cancel. Higher power terms here thus do not need any dimensionful prefactor. 
Moreover, there are no components of the inverse metric (which would be l/p\v) for our 
diagonal triads) in contrast to what is required in higher curvature terms. 



Curvature expansion. To see how this is reconciled, we expand the Hamiltonian explic- 
itly in io after writing kj = ioki- This corresponds to a slowly varying field approximation 
with respect to the lattice size. For the (+/, + J)-plaquette, a single term in the sum fl2T|) 
becomes 

'^\i,j^v,jCv+i,j + 2s^j jCy^js^+i^j = j^ilkjkj + ^l^^l (kikj,j + kjkjj + 2kjkij^ 

kjkjjj + 4{kikjjj + kjkjjj + kjjkjj + kj^jkjj) 
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Figure 3: Four plaquettes adjacent to vertex v in the (/, J)-plane. The arrows indicate the 
directions in which the relevant holonomies are traversed. 



+2kjjkj^j - -i^kjkjik] + k]] 



f28) 



(Commas on the classical field kj indicate partial derivatives along a direction given by 
the following index.) For a fixed direction K there are in total eight terms to be included 
in the sum ( 12T|) . They are obtained from ( l28l) by taking into account the four plaquettes 
in the (J, J)-plane meeting at vertex v (Fig. [3]) and considering both orientations in which 
each plaquette can be traversed. While the latter merely boils down to symmetrization 
over I and J, the former requires some care, noting that in the Hamiltonian constraint 
(I2T1) h^-i means h~^j j. The contribution (l28l) corresponds to plaquette 1 of Fig. [3] and 
has (Ji = aj = Ok = 1. Accounting for the overall sign dictated by the cr-factors, one can 
obtain the expressions for the three remaining plaquettes 2, 3 and 4 following the recipe 
below. 



plaquette 


extrinsic curvature components 


sign 


(1) 


ki{v + I/2) 


kj{v + I + J/2) -ki{y + I/2 + J) 


-kj{v - 


VJ/2) 




(2) 


~kj{v - 1/2) 


kj{v-I+J/2) ki{v-I/2 + J) 


-kj{v - 


VJ/2) 


(-1) 


(3) 


-ki{v - 1/2) 


-kj{v -I- J/2) ki\v - 1/2 - J) 


kj{v - 


J/2) 




(4) 


ki{v + I/2) 


-kj{v + / - J/2) -kj{v + 1/2 - J) 


kj{v - 


J/2) 


(-1) 



The first column designates a plaquette number, whereas the last one indicates the overall 
sign factor. The other four columns show the correspondence between the relevant link 
labels. 

After the symmetrization over all four plaquettes (traversed in both directions), the 
cubic terms drop out 



I'iik.kj 



^-^kjkj{k] + fc}) 



,,2/^4 



(29) 



+ ^ {kikj^jj + kjki^n + 2(^/fcj,// + kjki^jj + kj^ikj^i + ~ki,jkj,j)^ + 0(£^) . 



Note that the link labels k were introduced as values of the extrinsic curvature components 
evaluated at midpoints of edges in the continuum approximation (!23|) of our basic non-local 
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variables. The expression above is written in terms of just two components kj{v) and kj{v) 
(and their partial spatial derivatives) Taylor expanded around the vertex v. The first term, 
when combined with B^ k and summed over all triples UK, reproduces the correct classical 
limit of the constraint Hk- This limit is obtained in two steps: we first performed the 
continuum approximation by replacing holonomies with mid-point evaluations of extrinsic 
curvature components. This would still give us a non-local Hamiltonian since each vertex 
contribution now refers to evaluations of the classical field at different points. In a second 
step we then Taylor expanded these evaluations around the central vertex v, which gives 
a local result and corresponds to a further, slowly-varying field approximation. 

Comparison with higher curvature terms. Here, the factor in the leading term 
together with a factor from B^^k through ( |T^ combines to give the Riemann measure 
of the classical integral. Higher order terms, however, come with additional factors of £o in 
(!29|) which are not absorbed in this way. The result is certainly independent of coordinates 
since the whole construction (l26|l in terms of kj is coordinate independent. But for a 
comparison with higher curvature terms we have to formulate corrections in terms of kj 
and as these are the components of classical extrinsic curvature and densitized triad 
tensors. Higher order terms in the expansion are already formulated with kj in coordinate 
independent combinations with ^o-factors. It remains to interpret the additional io factors 
appropriately for a comparison with low energy effective actions. 

This can be done quite simply in a way which removes the above potential discrepancies 
between our expansions and higher curvature terms in low energy effective actions. We 
simply use ([6]) to write £q = /p^ which is the only well-defined possibility to express £o 
in terms of the fields. Thus, inverse metric components 1/p^ directly occur in combination 
with kj factors as required for higher curvature terms. The fact that the cubic term in 
io in fl29|) drops out is also in agreement with higher curvature corrections since in that 
case only even powers of the length scale ip occur. Moreover, there are now factors of p^ 
multiplying the corrections. These are basic variables of the quantum theory determining 
the fundamental discreteness. Thus, factors of the Planck length occurring in low energy 
effective actions are replaced by the state specific quantities p^ . While the Planck length 
ip = a/ Gh is expected to appear for dimensional reasons without bringing in information 
about quantum gravity (it can just be computed using classical gravity for G and quantum 
mechanics for h), the p^ are determined by a state of quantum gravity. If expressed through 
labels fiyj, the Planck length also appears, but it can be enlarged when fi^ j > 1. Moreover, 
the lattice labels are dynamical (and in general inhomogeneous) and can thus change in 
time in contrast to ip. Although the form of corrections is analogous to those for low energy 
effective actions, the conceptual as well as dynamical appearance of correction terms is thus 
quite different. 

The terms considered so far could not give rise to higher time derivatives of the spatial 
metric. In general, higher time derivatives describe the effect of quantum variables ([3]) of 
the field theory, which appear in the expectation value of the Hamiltonian constraint in 
a generic state. Quantum variables are thus, in a certain sense, analogous to higher time 



16 



derivatives in effective actions which indicates that the correction terms they imply 
should combine with those in (l29l) obtained by expanding sines and cosines to higher 
powers of space-time curvature components. All corrections of these types should thus 
be considered together since they will eventually be mixed up despite of their different 
derivations. A computation of terms containing quantum variables requires more detailed 
information about the expectation value of the constraint operator in an arbitrary state. 
These terms are thus more difficult to compute, which also makes an interpretation of the 
remaining higher curvature terms alone, especially concerning their possible covariance, 
more difficulto We will thus focus from now on on corrections coming from commutators 
Bv,K to quantize inverse powers which are independent of the higher order corrections and 
even give rise to non-perturbative terms. Moreover, in Sec. l5.2.2] below we will demonstrate 
that those corrections are expected to be dominant in cosmological perturbation theory. 

3.1.3 Inverse triad corrections 

A direct calculation using ([8]) and (fT2!) shows that B^^k commutes with all flux operators 
and thus has flux eigenstates as eigenbasis, as it happens also in homogeneous models 
The action 

^ 1 /2 / 

B^^K\---,fiv,K,---) ■= (27r7£p) y\f^v,i + IJ'v,-i\\l^v,i + IJ'v-j\ (30) 



X \fJ'v,K + l^v-K + 1| ~ Y \l^v,K + fJ'v,-K ~ ^\ j I • • • ' ^^'u,K, ■ ■ •) 

directly shows the eigenvalues which do not agree exactly with the classical expectation 
exiv) = ^/\p^ {v)p-^ (v) / (v) I ~ \/\f^v,ifJ'v,j/ I-J'v,k\ (indices such that e^'-'^ = 1) for the 
co-triad (|T3l) which appears as a factor in the Hamiltonian constraint. But for large values 
^ 1 the classical expectation is approached as an expansion of the eigenvalues shows. 
Inverse triad corrections are obtained by extracting the corrections which B^^k receives 
on smaller scales. We introduce the correction function as a factor a^^K, depending on 
the lattice labels fi^j, such that By x = C(v,k€'v,k and ay^K ^ 1 classically, i.e. for fi^^^ ^ 
1. Comparing the eigenvalues of B^^k with those of flux operators in the combination 
\/\J^v,i^v,j/J^v,k\^ we find 



^It is sometimes tempting to "sum the whole perturbations series" of higher order terms simply by 
using the left hand side of (pS)) directly in effective equations without an expansion. However, this is in 
general not a consistent approximation since arbitrarily small higher order terms are included but other 
types of corrections such as higher time derivatives or quantum variables are completely ignored. There 
is currently only one known model where the procedure is correct since all quantum variables have been 
shown to decouple from the expectation values [20]. But this model, a free isotropic scalar in a certain 
factor ordering of the constraint operator, is very special. Under modifications such as a scalar potential 
quantum variables do not decouple and their motion implies further correction terms in effective equations 
not captured in the trigonometric functions arising from holonomies. 
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After having computed the operators and their eigenvalues, we can speciahze the cor- 
rection function to perturbations of the scalar mode. We reduce the number of independent 
labels by imposing fi^j + fj^v-i = l^v,j + fJ^v-j for arbitrary / and J. This corresponds to 
a metric proportional to the identity Sab for a scalar perturbation. We then assign a new 
variable p{v) = 27i'yip{fiyj + fiv,-i) to each vertex v, which is independent of the direction 
of the edge / and describes the diagonal part of the triad. Quantum numbers in eigenvalues 
of the lattice operators can then be replaced by p{v), and the resulting functions compared 
with the classical ones. The remaining subscript v indicates that the physical quantities are 
vertex-dependent, i.e. inhomogeneous. Then the averaging over the plaquette orientations 
in the constraint becomes trivial and the total correction reads 

a, = (Vl/i. + 1/2| - Vl/i. - 1/2|) (32) 

i.e. 

= (^^\p{v) + 27r^il\ - ^\p{v)-2n^il\^ . (33) 

We will continue analyzing these correction functions in Sec. 13.31 after having discussed 
how such functions also enter the spin connection and matter terms. 



3.1.4 Spin connection 

So far, the holonomies we used only contributed the extrinsic curvature terms to the 
Hamiltonian but no spin connection terms at all. In the procedure followed here, we thus 
have to quantize r^[£'] directly which is possible in the perturbative regime where line 
integrals of the spin connection have covariant meaning. This gives rise to one further 
correction function in the effective expression of the spin connection 

= -e^^'ej (^a[,e5 + ^eleid^Ai^ , (34) 

as it also contains a co-triad ( |T3l) . Since the triad and its inverse have a diagonal form 



\/\detE\ 
with the components given by 



e, = e(,)<5l (35) 



= —£== = {eiy\ detE=pyp^, (36) 
det E\ 

the spin-connection simplifies to 

= ere(^)9,e(,). (37) 
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In terms of components of a densitized triad it reads 



Since there are many alternative choices in performing the quantization of such an ob- 
ject, but not much guidance from a potential operator in the full theory, we first discuss 
general aspects one can expect for the quantization of the spin connection in a simple ver- 
sion. It includes corrections of inverse densitized triad components by correction functions 
in each term of (l38l) . We thus mimic a quantization to the extent that expectation values 
of classical expressions containing inverse powers of p acquire a correction factor 

' '*(") (39) 



p^{v) p^{v) ' 

where the correction functions (3j are kept different from the function a used before because 
the object to be quantized is different. There will also be corrections from the discretization 
A I of partial derivatives iodj, but we ignore them in what follows for the same reason which 
allowed us to ignore such effects from the loop holonomy quantizing F^^. The effective 
analog of (138|) is then of the form 

At this stage the triad components, corresponding to different orientations, can be put 
equal to each other in effective equations, p^ = p'^ = p^ = p. This implies an analogous 
relation between the correction functions Pi = Pj = Pk = Po- Comparing (HO!) with 
the ansatz T} = P^, we conclude that also the spin connection receives a correction 
function P = Pq. 

For a precise quantization we observe that we need terms of the form ^qFJ^F^ and 
i^daTl in the constraint since one factor of the Riemann measure will be absorbed in 
the commutator Byj. To quantize ioT^^, we combine io with the partial derivative dj in 
fl37|) to approximate a lattice difference operator Aj defined by (A//)^ = f^+j — for 
any lattice function /. A well-defined lattice operator thus results once a prescription for 
quantizing the inverse triad has been chosen. One can again make use of Poisson identities 
for the classical inverse which, however, allows more freedom than for the combination of 
triad components we saw in the Hamiltonian constraint. Such a freedom, corresponding to 
quantization ambiguities, will also be encountered when we consider matter Hamiltonians. 
For any choice we obtain a well-defined operator which would not be available without the 
perturbative treatment since the full spin connection is not a tensorial object. 

An explicit example can most easily be derived by writing the spin connection integrated 
along a link e^j as it might appear in a holonomy, 

Je,,i V|detE| 
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using the lattice difference operator A/ ^ iodj. We then have to deal with the inverse 
powers explicit in the fraction and implicit in the co-triad ej. The latter is standard, 
replacing e/ by io^hi{hJ^ ,Vy} based on (fT3l) . The inverse determinant in the fraction 
cannot be absorbed in the resulting Poisson bracket because (i) it does not commute 
with the derivative and (ii) absorbing a single inverse in a single co-triad would lead to a 
logarithm of in the Poisson bracket which would not be well-defined. It can, however, be 
absorbed in the flux i^p^ if we do not use the basic flux operator F^ k but the classically 
equivalent expression 

Fy,K ^ V = \(^l5tKfH,e^''eyj 

= -l(47r7G)-2^ ^i^J^''''HnKMhj\Vy}hj{h-j\Vy}) (41) 

IJ <T/e{±i} 

which is analogous to expressions used in [2T] . Since there are two Poisson brackets, we can 
split the inverse evenly among them, giving rise to square roots of K, in the brackets: 



^y\detE\ K 



J2 Yl ^i^je''''tT{T^K)hi{hj\y^,}hj{h-j\^/V^}){A2) 



' IJ CT/e{±i} 

The remaining factor of io is absorbed in e/ inside the derivative which is quantized follow- 
ing the standard procedure. A well-deflned quantization of spin connection components 
thus follows, which is not local in a vertex since the difference operator connects to the 
next vertex. Similarly, the derivative of the spin connection needed in the Hamiltonian 
constraint leads to further connections to next-to-next neighbors. 

Explicitly, one can thus write an integrated spin connection operator quantizing ^ := 
dte?r^ as 

Replacing the commutators by classical expressions times correction functions a deflned 
as before and ai/2 deflned similarly for a commutator containing the square root of the 
volume operator leads to an expression 

(r})efT = «i/2(p>i/2(/)e/^'e(^)9e(«(/)e,) 

= «i/2(p>i/2(p')«(p')ri + «i/2(p')«i/2(/)«'(p')e/e/'^e('^)a,p(^) 
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where the prime denotes a derivative by . Using the relation p'^ = eiex whenever 
ejiK = 1 between densitized triad and co-triad components allows us to write 

(r})efr = ai/2{p')ai/2{p^)a{p^)ri + ai/2{p')ai/2{p^)a\p^)eie/''e^''^dc{ej^ 

= «i/2(p>i/2(/) ( a(/)ri + J2 c^'ipVrU (44) 

for the effective spin connection components. For scalar modes, using that all p^ at a given 
point are equal, this can be written with a single correction function 

f3[p{v)] = a^/2[p{v)]\a[p{v)] + 2pa'[p{v)]) (45) 

for F}, where a' = da /dp. 



3.2 Matter Hamiltonian 

Matter fields are quantized by similar means in a loop quantization, using lattice states, and 
then coupled dynamically to geometry by adding the matter Hamiltonian to the constraint. 
For a scalar field 0, the momentum vr = a/| det E\^/N is a density of weight one. In the 
0-representation, states will simply be of the form already used for the gravitational field, 
except that each vertex now also carries a label i/^, G M describing the dependence on the 
scalar field 0(f) through exp(iz/^0(f )) [22]. Well-defined lattice operators are then given 

by exp(ii/o0t;); ^oi any z/q G M, which shifts the label by z/q. The momentum, with its 
density weight, has to be integrated before it can meaningfully be quantized. We introduce 

:= / d^XTX ^ (-I'niy) 

J Ry 

where is a cubic region around the vertex v of the size of a single lattice site. Since we 
have {(j){v),Pw} = Xn^i'^) terms of the characteristic function Xfi(^) = 1 if f G i? and 
zero otherwise, a momentum operator P„ must have eigenvalue hu^ in a state introduced 
above. 



3.2.1 Inverse triad corrections 



For the matter Hamiltonian of a scalar field with momentum vr and potential f/(0) we 
have the classical expression 



d^xN{x) 



2,/d^ 



n{xf + 



2^/detq 



da(l){x)db(l){x) + v/d^C/(0) 



containing inverse powers of the metric, too. It can be quantized by loop techniques 
making use of identities similar to f|T3l) . One first generalizes the identity to arbitrary 
positive powers of the volume in a Poisson bracket. 



Ka,y;} = 47r7GrK:-V 



(46) 
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and then combines such factors with suitable exponents r to produce a given product of 
triad and co-triad components. Since such identities would be used only when inverse 
components of densitized triads are involved and a positive power of volume must result 
in the Poisson bracket, the allowed range for r is < r < 2. Any such Poisson bracket will 
be quantized to 

using holonomies hyj in direction I with tangent vector e^. Since holonomies in our lattice 
states have internal directions Tk for direction K, we can compute the trace and obtain 

KR. = ^^^atrir^KMKiK, = ^^{B^}, - B^k^k, (47) 

where, for symmetry, we use both edges touching the vertex v along direction K and B^j j^ 
is the generalized version of (1251) : 

■■= {s.,kV;c.,k - c.,kV:s.,k) (48) 

The exponent used for the gravitational part was r = 1, and r = 1/2 already occurred in 
the spin connection, while the scalar Hamiltonians introduced in [23112^. which we closely 
follow in the construction of the matter Hamiltonian here, use r = 1/2 for the kinetic term 
and r = 3/4 for the gradient term. With 

^abc^ f Ai T/l/2ir/ij T/l/2ir/ifc t/1/2i fO-rr^,n\3 ^abc^ ^a^b^c Q{27r'yG)^ 

e eijk{A^, 1/ ' }{Al, VJ VJ \ = (27r7G) e eijfc— ^ = ./^ 

for a lattice site volume ^ i\\ det(e^)| and 



one can replace the inverse powers in the scalar Hamiltonian as follows: For the kinetic 
term, we discretize 

r2 „(n,\'2. r> 2 



Then, the classically singular 

^ = ( = ( t^tSt^.^^'-^^AAI K'/'HAl K'/'HA',, V^''} ) (49) 



will be quantized to 



f 1 |'r(1/2) R(l/2)^^^ { uC^/"^) R(l/2)^^i / r(1/2) R^^/^) n efc \ 



2 
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Similarly, we discretize the gradient term by 

J Vdet^ ^ Vdet q{v) ^ K 

where we replace spatial derivatives da by lattice differences A/. Now, using 

we can quantize the metric contributions to the gradient term by 



1 

.JMN^ ^r(3/4) r(3/4) X ^6(3/4) r(3/4)n™ 



JKL^ |'r(3/4) -(3/4) s.j /A(3/4) R(3/4)x.fc 



In addition to the fact that we are using different values for r in each term in the grav- 
itational and matter parts, giving rise to different correction functions, the matter terms 
are less unique than the gravitational term and can be written with different parameters 
r. This corresponds to quantization ambiguities which will appear also in effective equa- 
tions and which could have phenomenological implications. Some choices are preferred 
since they give rise to simpler expressions, but this does not suffice to determine a unique 
quantization. Instead of using r = 1/2 in the kinetic term, for instance, we can use the 
class of relations 



(dete)'^ /I 



-e 



./\d^\ (detE)('^+i)/2 V6 

X{A* , •l/(2'=-l)/3*=}{^J^ y{2fc-l)/3fc||^i^ y(2fc-l)/3fc| 

for any positive integer k to write the inverse determinant through Poisson brackets not 
involving the inverse volume (see also the appendix of [25]). This determines an integer 
family of quantizations with = {2k — l)/3k > |. For = 2 we obtain the previous 
expression, but other choices are possible. Moreover, using the same r in all terms arising 
in gravitational and matter Hamiltonians can only be done in highly contrived ways, if 
at all. There is thus no clearly distinguished value. From now on we will work with the 
choices specified above. 

On regular lattice states, all ingredients are composed to a Hamiltonian operator 

HM = E'V. E -i-2-3."'^Bri<'5C':i-) (51) 

IJK,aie{±l} I 
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Figure 4: Behavior of the correction function a. It approaches one from above for large 
arguments. For small arguments, the function is increasing from zero and reaches a peak 
value larger than one. Also shown is the limiting case r = 2 which does not show a peak 
but a constant correction function for > 1. 



3.2.2 Matter correction functions 

As before, we compute eigenvalues of the operators 



5i J := (27174) ^ ^^'^"'"'-^ (52) 

where the subscript of the volume operator indicates that its eigenvalue in a lattice state 
is computed according to (IT^ with a shifted label of link e^_/. Thus, the eigenvalues are 

(53) 

compared to the classical expectation 

(27774) ^"V*;,/ + ^v-l\'''^\[i-v,J + [i-v-j\'''^\[i-v,K + ^v-kX'"^'^ 

for V^'^tK- For any r, correction functions 



result. The main examples of r are seen in Fig. HI 

Imposing again that ii^j + liv-i = fJ'v,j + l^v-.j for arbitrary I and J and introducing 
piy) = 2'7i'yip{fiyj + fiv-i), we obtain the effective correction functions 

a^'-^piv)] = b(t;)r-'-/2(b(t;) + 27774 r/'-b(^)- 27774 r^') ■ 
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This can be used to write the effective matter Hamiltonian on a conformally ffat space 
Qab = \pix)\Sab as 

where comparison with fl5T|) shows that 

D[p{v)\ = a^^/^\p{v)f and a[p{v)] = a^'"^\p{v)f . (55) 



3.3 Properties of correction functions from inverse powers 

We have derived several different correction functions, making use of different parameters 
r. In most cases one could make different choices of such parameters and still write the 
classically intended expression in an equivalent way. This gives rise to quantization ambi- 
guities since the eigenvalues of 5^^]^ depend on the value r, and so will correction functions. 
In addition to the ambiguities in the exponents r, one could use different representations 
for holonomies before taking the trace rather than only the fundamental representation 
understood above [26l|27j. In this case, we have more generally 



Eigenvalues of such operators can be expressed as 
3 

B^vK = ■/ ■ , ■ I i\ (27174) \y^v,i^y^v-i\'^'^\^v,j^y^v-j\'''^ m\y,^^K^^„_K^2m\'l'^ 
rj{j + l)(2j + 1) ^ J^. 

(56) 

which leads to the general class of correction functions 

3 r ^ 

(^v,K = —T—^r777--^\f^v,K + V m \fi^,K + f^v,-K + 2mf/^ . (57) 

rj(j + l)(2j + 1) J^. 

After imposing isotropy the last expression becomes 



a 



rj(j + l)(2j + 



— lyul"^ 2 m|yU + m|^''^ (5J 



m=-j 



which is shown for a few cases in Fig. [5] 

For large j, the sum in a^"^'^^ can be approximated using calculations as in [27j. The idea 
is to consider two cases: i) /i > j and ii) fi < j separately. In the former, the absolute values 
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Figure 5: Behavior of the correction function a for larger j. The general trend is similar 
to the case for j = 1/2, but there are [j] + 1 spikes at yU = 1, . . . , j for integer j and 
[1 = 1/2, ... ,j otherwise. Above the peak, the function is smooth. 



can be omitted as all the expressions under the sum are positive. Then the summation is 
to be replaced by integration to yield 



a 



3,-,l- 



rj{j + l){2j + l) 



r + 4 
~r + 2 



((^+l)i+^-(/2-l)i+^) 
((/i + l)i+^-(/i-l)i+^) 



/i > 1 



where fi := fi/j. In the second case, the terms in the sum corresponding to m < fi and 
m > n should again be considered separately. The end result, however, is very similar to 
the previous one 



a 



rj(j + l)(2j + l) 



1 



r + 4 
r + 2 



((/2 + l)i+^-(l-A)^+^) 
((/i+l)i+^ + (l-^)i+i) 



, fl<l 



After some rearrangements and using that j ^ 1 these two expressions can be combined 
into a single one as 



a 



6/2 



r{r + 2)(r + 4) 



+ l)i+i(r + 2 - 2/2) + sgn(/2 - l)|/i - l\i+\r + 2 + 2/i)] .(59) 



The approximation is compared to the exact expression of the correction function ob- 
tained through eigenvalues in Fig. [61 As one can see, the spikes are smeared out by 
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Figure 6: Comparison between the correction function (1571) and its approximation (l59l) . 
The spikes are smeared out by the approximation. 



the approximation (except for the point fx = 1 where the approximation remains non- 
differentiable at second order which is not visible from the plot). The general trend, how- 
ever, is reproduced well even below the peak. For applications in effective equations we 
note that the approximation might be considered more realistic than the exact eigenvalue 
expression because those equations would be based on semiclassical states. Since such 
states cannot be eigenstates of the triad but must only be peaked on a certain expectation 
value, they will automatically give rise to a smearing-out of the spikes in the eigenvalues 
as discussed in more detail in Sec. I4.1[ 



3.3.1 Asymptotic behavior 

This class of correction functions parameterized by two ambiguity parameters r and j 
captures the most important general properties of such functions, including the position of 
their maxima at /i ~ 1 (or /i ^ j) and the initial power law increase for small /i (determined 
by r) [271 [28] • It is indeed easy to see that all correction functions have the correct classical 
limit on large scales, such as 

« 1 , _L_ (r-2)(.-4)4(3,^ + 3,-l) ^ ^ ^ 

for fl32|) . Moreover, for small the correction function goes to zero as 

a(^'^)(/i) ^ (2/i)2-i, (61) 

which ensures boundedness of the quantized co-triad e(^K) oc a^/Ji oc fP' (using r = 1 for 
this case as in f[T51) ). when /i ^ 0. The same is true for higher j since evaluated at = 
the sum of odd terms gives zero. 
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This function is not smooth but has a cusp at its maximum at /i = 1/2, or more 
generally a cusp at integer or half- integer values between and j. The second derivative 
a" is always positive while a' changes sign between any two cusps. To the right of the cusp 
at the largest /i, the derivatives satisfy 

a' < 0, a" > . (62) 

Note that the approximation used for larger j smears out the cusps and does not everywhere 
have positive second derivative. The behavior above the peak and the general increase 
below is, however, reproduced well by the approximation. The definite sign of a" has far- 
reaching implications in the quantum corrected equations of motion [6]. Small corrections 
then add up during long cosmic evolution times which would not be realized if, e.g., a 
would oscillate around the classical value which is also conceivable a priori. 



3.3.2 Small-scale behavior and ambiguities 

We will mostly use here and in cosmological applications of the corrected perturbation 
equations of [5J the behavior for larger values of /i above the peak. On very small scales, the 
approach to zero at /i = is special to operators with U(l)-holonomies as they appear in the 
perturbative treatment here. In particular, as we have seen explicitly the volume operator 
V and gauge covariant combinations of commutators such as tr(r*/i[/i^^, V]) commute. It is 
thus meaningful to speak of the (eigen-)value of inverse volume on zero volume eigenstates. 
For non-Abelian holonomies such as those for SU(2) in the full theory, the operators become 
non-commuting [29]. The inverse volume at zero volume eigenstates thus becomes unsharp 
and one can at most make statements about expectation values rather than eigenvalues 
which again requires more information on semiclassical states. Then, the expectation values 
are not expected to become sharply zero at zero volume, as calculations indeed show [30] . 
In addition, also here quantization ambiguities matter: We can write volume itself, and 
not just inverse volume, through Poisson brackets such as [29] 

After a lattice regularization as before, using \4 ^ £q| det(ejj)|, we obtain 

2 / i j k ^ 



( abc 

whose quantization, making use of commutators, differs from the original volume operator 
( IT2|) . If non-Abelian holonomies are used, it would not commute with the full volume 
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operator of [31] or [32]. This clearly shows that the usual quantization ambiguity also 
applies to what is considered the relevant geometrical volume. (Related ambiguities for 
flux operators have been discussed in [21].) It is not necessarily the original volume op- 
erator constructed directly from fluxes, but could be any operator having volume as the 
classical limit. For finding zero volume states to be related to classical singularities, for in- 
stance, dynamics indicates that volume constructed in the more complicated way through 
commutators with the original volume operator is more relevant than the volume operator 
constructed directly from fluxes [29]. Thus, specific volume eigenstates have to be used 
with great care in applications with non-Abelian holonomies. Also, the behavior of correc- 
tion functions below the peak value, especially whether or not they approach zero at zero 
volume, is thus less clear in a general context. In any case, below the peak positions scales 
are usually so small, unless one uses larger j, that perturbation theory breaks down. The 
behavior above the peak, by contrast, is robust and gives characteristic modifications to 
the cosmological evolution of structure. 

4 Effective Hamiltonian 

Calculations of distinct terms in the constraint presented in the preceding section can now 
be used to derive effective Hamiltonian constraints. 

4.1 Expectation values in semiclassical states and quantum vari- 
ables 

The derivation of an effective Hamiltonian constraint proceeds by computing expectation 
values of the constraint operator in semiclassical states which are superpositions of our lat- 
tice states peaked on perturbative metric and extrinsic curvature components. Such states 
are easily constructible although, for the order we are working at here, we do not need to do 
so explicitly. The peak values of perturbative fields are thus in particular diagonal which 
means that expectation values can easily be computed via Abelian calculations]^ The only 
complication arises from the fact that we are necessarily dealing with operators as products 
of holonomies and fluxes which are not simultaneously diagonalizable. It is most conve- 
nient to use the triad eigenbasis | . . . ,Hv,i, ■ ■ •) for triad or inverse triad operators, and a 
holonomy eigenbasis for products of holonomies. This was implicitly assumed previously in 
the curvature expansion and when using inverse triad eigenvalues for correction functions. 
However, for expectation values of the complete constraint operator as a product of holon- 
omy and co-triad terms we need to transform between the two eigenbases, which as usually 

^Although initially SU(2)-holonomies appear in the constraint, they only refer to lattice directions such 
that those holonomies are of the form exp(^r/). While these matrices do not commute among themselves 
for different /, one can easily re-arrange the order; see, e.g., [33l[3] for a discussion of the analogous effect 
in symmetric models. The special form of SU(2)-matrices occurring in our context is also the reason why 
we can take the trace in the Hamiltonian constraint explicitly and reduce it to a product of sines and 
cosines of curvature components. 
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is possible by inserting sums over complete sets of states: {iplHltp) = J2i{\'^Hi\I){I\H2\i') 
if {]/)} is the complete set of states and H = H1H2 is factorized into the two parts men- 
tioned above. For a complete treatment we thus need to compute matrix elements of Hi 
and i^2, not just eigenvalues. 

Nevertheless, the calculations presented before already provide the main terms under 
the following approximation: We assume, without loss of generality, that the complete set 
of states {]/)} contains a state we are interested in. More crucially, we assume that 
the spread a oi ip m basic variables is small. Under this assumption, i = 1,2, 

are dominated by (tplHilip) since (i) there is not much overlap with most other states in 
the complete set and (ii) the states {ip), having small spreads, are as close as possible 
to eigenstates of Hi and H2, respectively. With Hi being a product of holonomies and 
H2 depending on fluxes, the spreads required in this construction cannot be arbitrarily 
small because they are restricted by uncertainty relations. This implies that additional 
corrections not computed before arise due to the unavoidable spread of semiclassical states. 
As a direct consequence of spreading, such terms depend on parameters such as a which 
are nothing but the quantum variables (I3l) mentioned before. These variables necessarily 
feature in a complete effective Hamiltonian, describing how spreading and deformations of 
the state back-react on the peak position [T3]. 

Apart from these quantum variable terms, the main effective Hamiltonian then is of 
the form (tplHilip) {■ip\H2\ip) where is a semiclassical state peaked on a given classical 
geometry. Each of the two factors is of the form 

5^0^|^(^")p~5^0^/(/x-/xo)~ / d/iO(/x)/(/x-/io) = (0^/)M (63) 

where are eigenvalues of an operator O and is the peak value of the state \ip) in the 
basic variable /i whose eigenbasis is used. On the right hand side, we see that the effect of 
computing an expectation value in a semiclassical state is mainly, to the given order, that 
eigenvalues appear in a form convoluted with the shape of the semiclassical state. 

In such a convolution, sharp features in eigenvalue functions such as the spikes in Fig. |5] 
will be smeared out. But otherwise the general behavior is already displayed well by 
explicit eigenvalues, and, similarly, higher order curvature corrections are close to what we 
computed before. For general features we can thus avoid dealing with details of states and 
their convolution with eigenvalues. The form of effective Hamiltonians we arrive at is thus 

Hes[N] = ^J^d'xN{x)a[p](^~3k^ + AK (64) 
+ j &'xN{x)Aq 
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for metrics including scalar perturbations in longitudinal gauge. Note that the correc- 
tion functions depend functionally on the field p{x), not p, which shows that their scale 
is uniquely determined by the state irrespective of choices of coordinates. Unspecified 
correction terms are higher order curvature corrections A;^ (see Eq. (1291) ). discretization 
corrections A^^'^^'' from different spatial derivative terms in the constraint, and terms con- 
taining quantum variables Aq which arise from metric as well as matter fields. 

This form of effective constraints also demonstrates potential effects of using SU(2) 
representations different from the fundamental one. Notice that we did not compute this for 
the higher curvature expansion since the required traces of different Pauli matrices are more 
involved. But it is clear that this can only change the coefficients in the expansion A^^- since 
it always remains at a perturbative level. Generally, larger values of j mean that curvature 
corrections will become important at smaller curvatures compared to j = 1/2, and thus 
coefficients in an expansion will increase with j. The effect in inverse triad corrections a[p], 
which we did compute explicitly here, is more pronounced since j determines the scale at 
which one enters the non-perturbative regime of such inverse triad corrections. The main 
difference between larger values of j and the minimal one is that in the former case peaked 
states exist whose spread is smaller than the peak position of eigenvalues of an inverse 
triad operator. When this is realized, the non-perturbative branch of increasing behavior 
between /i = and /i = j is not completely washed out in the convolution but remains 
visible in an effective Hamiltonian. 

In an effective Hamiltonian this consequence is obvious, but it was not always clear 
from the underlying difference equation in isotropic models. There, the discrete stepsize, 
determined by the SU(2)-representation of holonomies in the constraint, equals the peak 
position of inverse triad eigenvalues (see, e.g., [M]). Thus, the discrete step in the difference 
equation always jumps from zero directly to the peak when the same representation is used 
for all holonomies occurring in the gravitational as well as matter parts of the constraint. 
One could thus argue that dynamics will be insensitive to the value of the representation. 
Effective equations, if they are applicable in this small-scale regime, show that this is not 
necessarily so. The representation enters higher curvature terms differently from inverse 
triad corrections, thus allowing effects of the non-perturbative part to remain. 

4.2 Technical issues 

We now illustrate some of the more important choices we made in constructing the con- 
straint operators used here. 

4.2.1 Quantization procedures 

Our construction is suitable for a treatment of cosmological perturbations within loop 
quantum gravity, but it does circumvent some of its general aspects. First, we do not take 
into account full non-abelian features; they can be included perturbatively but are not 
required for our selection of mode and gauge. Secondly, we do not allow irregular lattices 
or valence higher than six. Also this can be included by summing operator contributions 
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over triples of edges (as they are constructed in the full setting). Detailed coefficients in 
effective expressions may then change, but not qualitative effects. Moreover, as already 
mentioned the labels coming with additional edges or higher valent vertices are redundant 
for cosmological perturbations. 

We have presented higher power corrections using "holonomies" based on extrinsic 
curvature rather than connection components since this simplifies the calculations con- 
siderably. Using the background, it is mathematically possible to define such objects, 
although in a full background independent setting only holonomies of a connection would 
be well-defined. We use this mainly as a first possibility to demonstrate which types of 
corrections one expects and will discuss now how general the resulting expressions can be 
considered to be. This refers to corrections to terms of the Lorentzian constraint which, 
schematically, can be written as 

F + K^ = dA + A^ + {A-rf = d(r + K) + {r + Kf + (65) 

to be multiplied by triad components dealt with by Poisson brackets. 

Using extrinsic curvature as basic object, one obtains trigonometric functions of its 
components which when expanded give higher power corrections to dK + K^. But since 
the spin connection has been split off from the basic object, one has to quantize it indi- 
vidually and add suitable combinations for dF + to the constraint. Here, we assume 
that the cross-term TK does not contribute which is indeed the case for diagonal triads 
(implying anti-symmetric spin connection components) and extrinsic curvature. This is 
not much of a restriction: K is required to be diagonal for i^-holonomies to simplify the 
calculations. Moreover, the perturbative non-diagonal part of F must be antisymmetric 
because it perturbs an S0(3)-matrix. If there is a diagonal contribution, e.g. from a spa- 
tially curved background, it can be combined with K. As for corrections, we have higher 
power corrections in the quantization of dK + and inverse triad corrections in dF + F^ 
since the spin connection contains inverse triad components. 

Using A-holonomies gives, at first sight, a different picture. Now, F = dA + A^ receives 
higher power corrections, but the spin connection is not quantized directly, not giving 
immediate inverse triad corrections. One rather has to proceed as in the full theory [18\ 
where the term (A— F)^ = in the constraint is re-written using flTB]) . The double-Poisson 
bracket (fT6l) used to quantize extrinsic curvature now leads to additional corrections. In 
particular, since inverse triad quantizations have been used in H^^^ in (fT71) . corresponding 
corrections do arise which are qualitatively similar to those in a direct quantization of the 
spin connection. One thus expects similar types of corrections, as with i^-holonomies and 
F-quantizations, although in different combinations. 

In our construction, i^-holonomies arose from A-holonomies through a perturbative 
expansion in non-diagonal components. When using A-holonomies on a spatially curved 
background such as a closed Friedmann-Robertson- Walker model, it is furthermore neces- 
sary to take lattice effect of an inhomogeneous model [7] or related effects [35] into account. 
A fine enough lattice (A/" ^ 1) is required for a semiclassical expansion of holonomies since 

exp{iA) ~ exp(i(F + K)) ~ exp{i{K + Af~'^/^T + 6T)) 
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— 1 /3 ~ 

with the background spin connection possibly of the order T = Vq T = 0(1) can be 
expanded in all terms only if A/" is large. (For a closed isotropic model, for instance, 
r = 1/2 [36].) The number of vertices J\f enters through ioT = J\f^^^^T in holonomies. 
The spin connection perturbation 6r will be small in perturbative regimes such that we 
can always write 

exp{i{K + Af-^/^T + ST)) = exp{i{K + AT-^/^f ))(1 + iST + ■ ■ ■) . 

But the remaining exponential also has to reduce to the leading terms of an expansion in 
semi classical regimes. While K is then automatically small, this may not be the case for T. 
Without the reduction by J\f^^^'^ for a fine lattice, one could not expand the exponentials 
to reproduce the polynomials in K and F classically occurring in the constraints. 

We have focused here on the first type of quantization which is simpler to compute ex- 
plicitly but may not be as close to the full theory. Using A-holonomies, no diagonality can 
be used, but perturbative treatments of non-diagonal components are possible. It is thus 
possible, though more involved, to compute correction terms obtained through different 
quantization schemes and to compare their consequences, in particular those at the phe- 
nomenological level. A first result in this direction follows from the perturbation equations 
derived in [5] which show that effects of inverse triad corrections in F are less significant 
than those in the commutator [6]. Thus, one can hope that the precise quantization pro- 
cedure of curvature is not very important for physical aspects accessible so far. A detailed 
investigation of all consequences can nevertheless provide important guidance as to which 
procedure should be pursued in the full theory. 

4.2.2 Different types of correction functions 

In fact, we have included in the computation of perturbative effects four different correction 
functions a, (3, D and a. All of them come from inverse triad corrections. With all functions 
coming from the same type of modification, one may wonder why they should not all be 
identical. 

It is clear from the procedure that these functions arise from different classical func- 
tionals of densitized triad components. For instance, a comes from the antisymmetric 
part of EfE\/^\ detE\ while a comes from the symmetric part. They could be related to 
the same correction, but the quantization requires quite different rewritings (125!) and (!50|) 
of the corresponding terms in Hamiltonians such that correction functions will differ. In 
particular, they come with different parameters such as r. On top of that, each correction 
function is subject to quantization ambiguities. As we have seen, however, the typical 
behavior is robust under changes of the parameters. In particular, all correction functions 
have the same qualitative properties and differ only quantitatively in a way parameterized 
by a few parameters. 
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4.2.3 Implications for gauge issues 

The assumptions on states used to derive effective constraints liave a bearing on tlie gauge 
issue. By specifying the peak value of a spatial geometry and its extrinsic curvature in a 
semiclassical state we are fixing the spatial diffeomorphism constraint rather than solving 
it by averaging as done in the full theory [37]. Choosing the form of peak values partially 
implements a chosen gauge, but still allows some freedom. We also note that even though 
spatial diffeomorphisms are fixed, one still has to impose the constraint. This will give rise 
to one of the cosmological perturbation equations as is clear from [5]. 

Fixing the diffeomorphism constraint also implies a different viewpoint for the Hamilto- 
nian constraint operator of the loop quantization. In the full construction [18], one makes 
use of diffeomorphisms in order to make the operator more independent of the choice of 
edges used to quantize curvature. When diffeomorphisms are fixed, this is no longer pos- 
sible and effective constraints would depend on precisely how such edges are chosen. We 
have fixed this freedom here by laying the edges entirely on the lattice resulting in a graph 
preserving operator. Thus, holonomy corrections in the constraint depend on the lattice 
spacing provided through a state implementing the background geometry. While this sim- 
plifies the calculations without leading to significant quantitative changes in coefficients, 
we are as a consequence disregarding the creation of new vertices by the constraint opera- 
tor. Thus, N is constant for the construction, but may effectively be assumed to be slowly 
dependent on, e.g., the total volume (see [7] for more details). 

5 General implications for effective theory 

Our calculations, following the scheme to derive effective equations sketched before, have 
led to corrections (16^ which arise as leading order terms in an effective theory of pertur- 
bative loop quantum gravity. No complete expression has been derived, but characteristic 
terms are clear and lead to interesting phenomena [B]. Rather than studying one model 
in detail we have provided here an illustration of the general scheme: The characteristic 
feature of loop quantizations is the use of holonomies, which give rise to typical correction 
terms. They can be split into higher power corrections, which are always perturbative, 
and corrections to inverse powers of triad components which become non-perturbative at 
small scales. All these corrections are in addition to discretization and genuine quantum 
effects such as higher time derivatives. In this section we highlight conceptual conclusions 
that can be drawn from such a scheme. 

5.1 Basic variables in quantum gravity corrections 

Holonomy corrections arise through expectation values and thus depend on the basic vari- 
ables used in the quantization. Using commutators to quantize triad components, for 
instance, modifies the classical expressions in a way which can be computed through the 
explicitly available eigenvalues such as (jSSD of these operators. The occurrence of trigono- 
metric functions instead of direct curvature or connection components leads to higher power 
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terms when expanded in an effective constraint. Such corrections depend, by construction, 
on p{x)/ip = ilp{x)/ip and k{x) = iok{x), respectively, both of which are independent 
under rescahng the coordinates. They do depend, however, on the lattice size which de- 
termines the scales on which a state probes the field. 
If we split the fields into background parts 

p := — d^xp{x) and k .= — d^xk{x) (66) 
J yo J 

by integrating over a cube (sufficiently large to contain, say, the Hubble volume) of coor- 
dinate volume Vq = / d^x and perturbations 

6p{x) = p{x) — p and 6k{x) = k{x) — k (67) 



to set up cosmological perturbation theory [5] |j we can see from preceding constructions 
that it is not these fields directly which occur in correction functions. In isotropic loop 
quantum cosmology, the quantization is based on variables 



and 



k = V^"k - -57^ E ^oM^') = ^ E (69) 

which now appear as lattice averages in an inhomogeneous setting and provide the back- 
ground for cosmological perturbation theory. As such, they do not depend on the auxiliary 
coordinate volume Vq as they would in homogeneous models [17] but on the number Af of 
lattice vertices. These two quantities are related by Vq = Mil through the lattice size io, 
but M has significance as a parameter specifying the states rather than just being auxiliary 
as Vq. Similarly, basic variables of the inhomogeneous theory are function^ 

p{x) = ilpix) = —^p + 6p{x) = ^ E^(^) + ^P^^'^ ^'^^'> 
which directly occur in correction functions through fluxes, and 

k{x) = ioHx) = -^k + 6k{x) = ^ E '^(^) + '^'^(^) C^^) 
occurring in higher power corrections through holonomies. 



^In [Sj, only written here with a tilde have been used, but the tilde was dropped for notational conve- 



nience. 

7 



Note that io is used to rescale the inhomogeneous p{x) while Vq is used to rescale p as it is done in 
isotropic models [TT] . 
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This shows that the resulting equations are rescahng invariant when p, k and £o change 
simultaneously, a fact which was not always obvious in isotropic models based on the 
scale factor. As expected, the equations are also dependent on specifics (mainly Af) of 
the state whose dynamics is described effectively. This shows which states are suitable for 
perturbation theory and when perturbations break down. A perturbation scheme works 
only if 6p <^ p which from fITOl) implies that differences between local edge labels of the 
state (corresponding to 6p{x) = il6p{x)) must be small compared to the average lattice 
label J\f~^ p{v) (corresponding to the perturbative background value of p{x)). Since the 
labels are discrete, differences between them have a positive lower bound unless they are 
equal. Thus, the average label must be large compared to the discrete gap in the spectrum 
of labels. In our U(l)-theory, labels are integer valued which means that the average label 
must be larger than one, and local edge labels must not stray too much from the average. 
There is no such restriction from the curvature perturbations because curvature does not 
have a discrete spectrum. 

5.2 Quantum variables and classical limit 

Starting from the Hamiltonian (constraint) operator in any quantum theory, the quantum 
Hamiltonian is defined as a function on the projective Hilbert space determined by taking 
expectation values. This can be seen as the Hamiltonian function of a dynamical system 
whose phase space is obtained from the Hilbert space [HU [39l HO] . The system thus appears 
of classical form at least as far as dynamics is concerned, but each of its classical degrees of 
freedom is accompanied by infinitely many quantum variables ([3]). An effective description 
requires a further step, truncating the infinitely many quantum variables to a finite set 
[T2] . If this is done consistently, one obtains effective equations which amend the classical 
ones by quantum corrections. One often performs such a truncation by using a certain 
class of semiclassical states to compute expectation values of the Hamiltonian operator. 
The regime under consideration determines what a suitable set of semiclassical states is. 

Based on the assumed semiclassicality of states peaked at values p^ and kj, the expres- 
sions we derived give the main part of the effective Hamiltonian constraint computed as an 
expectation value in such states. Note that we did not explicitly compute expectation val- 
ues in states but read off corrections from operators by expanding trigonometric functions 
arising from holonomies or eigenvalues of inverse triad operators. Each of these correc- 
tions requires, strictly speaking, eigenstates of holonomies for higher order corrections in 
k or flux eigenstates for corrections as functions of p. But even if we were to compute 
expectation values in peaked states, the main corrections would be of the form read off 
from different eigenstates as seen by analogous calculations in the isotropic case pTl 1^ . 
In general, one has to use semiclassical states which are neither eigenstates of holonomies 
nor of flux operators. This gives rise to additional contributions depending, e.g., on the 
spread of the state. From the spread and other detailed properties of states one obtains 
contributions depending on additional independent quantum variables corresponding to 
fluctuations and correlations. For non-quadratic Hamiltonians or constraints, these quan- 
tum variables couple to classical variables and influence their motion. To some degree. 
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the appearance of additional independent quantum variables corresponds to higher deriva- 
tive terms in effective actions [13]. Thus, we obtain modified coefficients (from correction 
functions such as a), higher powers in momenta (from sinfc/ and cos/c/) and higher deriva- 
tive terms from quantum variables (interpreted as higher time derivatives) and from the 
discretization (higher spatial derivatives), which comprise all effects known and expected 
from effective actions. The first two arise as typical corrections by using holonomies. 

5.2.1 Basic variables vs. coarse graining 

In our treatment here we assumed that lattice scales are small compared to other scales 
of the relevant physical fields such as matter or classical metric modes to be obtained in a 
semiclassical limit. In such a context, it is sufficient to use the basic variables as they come 
as labels of a quantum state directly in effective correction functions. This is not possible 
in regimes where basic variables of the states are themselves strongly inhomogeneous as it 
necessarily happens when the discrete flux labels approach the lowest non-zero value 
one. Then, the perturbative condition Sfi^j <^ fi^j where 6fivj refers to the difference in 
nearby labels cannot be satisfied unless 6fiyj = 0, i.e. the labels are exactly homogeneous. 
Most likely, this happens in strong curvature regimes where perturbation theory would be 
expected to break down even classically. But since the discreteness of the labels plays a 
role in this simple argument, there can be regimes where classical perturbation theory would 
be applicable but the underlying lattice formulation would not seem to be in a perturbative 
regime. In such cases, one would have to coarse grain the basic variable, i.e. replace the 
basic lattice site variables by averages over larger patches of an intermediate scale. Then, 
the averaged labels would increase, relieving the contradiction between 6fiyj -C fiv,i and 
quantum discreteness. 

5.2.2 Orders of magnitude of corrections 

With several different correction terms, it is helpful to know whether in certain regimes 
some of them can be ignored. This can be difficult to determine in homogeneous models 
unless one makes special choices of ambiguity parameters such as large values of j [12] . In 
inhomogeneous situations it is often simpler to determine which corrections are expected to 
be dominant because they depend differently on the basic scale contained in p„ / ^ . These 
variables are parameters determining the state and thus the physical regime being probed. 
When p„ 7 is small, i.e. close to its minimum £p, inverse triad corrections are large. They 
decrease when p^j becomes larger, but this also implies larger and fewer lattice sites such 
that discretization effects become important. Moreover, in nearly isotropic configurations 
extrinsic curvature is given by 



as it follows from the Friedmann equation. The energy density scale p thus determines 
when curvature corrections are relevant. Since there is also a factor of p^j, curvature 
corrections increase with larger p^j just as discretization corrections. 




(72) 
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For a semiclassical regime we must have p^j > ip in order to reproduce closely the 
correct inverse powers of triad components. We must also have a discreteness scale pyj 
which is sufficiently small in order to avoid discretization effects already in, say, particle 
physics. This requires p^j to be much smaller than the typical physical scale squared, such 
as a wave length A of field modes or even the Hubble length a/ a. We thus have a range 
ip < pyj <^ or ip < pyj <^ (SvrGp)^^ if we express the Hubble length in terms of energy 
density. At the upper bound we ensure that discretization effects do not disrupt other 
physics used essentially in a given scenario. As a consequence of (1721) . this implies that 
higher order corrections in curvature are small, too. The dominant contributions are then 
given by inverse triad corrections which we have focused on in the preceding derivations. 
Note that the semiclassical range for pyj is large in late time cosmology, implying that 
corrections can be arbitarily small, for instance to the propagation of signals from gamma 
ray bursts. In the early universe, however, and in particular during inflation the energy 
scale is much higher, restricting the range more narrowly [B]. The best tests of quantum 
gravity effects are thus expected from early universe cosmology. 

5.2.3 Classical limit 

We have ignored in our calculations so far any detailed specifics of states and terms contain- 
ing quantum variables. Implicitly, we are thus assuming that such terms are subdominant, 
just as one assumes analogous terms to be subdominant in a derivative expansion of low 
energy effective actions. Under this assumption we reproduce classical expressions in the 
suitable limit, which proves that loop quantum gravity has the correct classical limit in 
this perturbative regime in the same sense as in usual effective theories. This statement 
certainly includes inhomogeneities in the perturbative sector considered here. For instance, 
the Newton potential and corrections on smaller scales can be obtained from perturbation 
equations derived from the effective constraints [6]. 

In effective theory, verifying the correct classical limit does not require one to construct 
explicit dynamical coherent states, not even approximately. This would certainly be of 
interest, but would be highly complicated and is rarely done in interacting field theories 
where one can nevertheless be certain about the correct classical limit. We emphasize 
that, in any case, a discussion of the semiclassical limit based on coherent states does 
require such states to be dynamical coherent states, or at least must involve statements 
on dynamical changes of state parameters to specify suitable regimes. This means that 
states must stay approximately coherent under evolution or, in a fully constrained theory 
such as gravity, solve the Hamiltonian constraint. If this is not realized, quantum variables 
and the back-reaction of spread and deformations on the classical variables are not under 
sufficient control to ensure the correct classical limit. 

There are two viable procedures to verify the correct classical limit of a quantum 
theory, be it a constrained or unconstrained system. First, one may use kinematical co- 
herent states to compute expectation values of the dynamical operators (a Hamiltonian or 
constraint operators) and then analyze the dependence of quantum variables in resulting 
equations of motion; if their effect on expectation values is small in suitable regimes, the 
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correct classical limit results. Secondly, dynamical coherent states can be used if they can 
be constructed at least approximately, which directly illustrates whether the dynamics of 
expectation values is close to the classical one. The second procedure is much more com- 
plicated for interacting theories since the full quantum dynamics would have to be solved 
at least approximately at the quantum level. The first procedure, by contrast, allows one 
to derive effective dynamical equations first and then approximate solutions to understand 
the behavior of expectation values. Thus, usually kinematical semiclassical states are used 
in explicit effective descriptions, followed by an analysis such as one in a derivative expan- 
sion in quantum field theory. Such a further analysis is always required when kinematical 
semiclassical states are used, and it can only be done in a regime dependent way to bring 
in conditions for when semiclassicality should be satisfied. We have done this implicitly in 
our discussion by assuming slowly varying fields as in usual derivative expansions, both in 
space by doing a continuum limit of the lattice states and in time by assuming quantum 
variables to be negligible. 

We emphasize again that even if one can demonstrate an "instant" classical limit by 
using kinematical coherent states, a dynamical statement would still require one to assume 
(or to show) that such back-reaction effects of quantum variables on expectation values are 
not strong. This picks the correct regime of states in which one has semiclassical behavior. 
Without such an additional analysis, kinematical coherent states would neglect the back- 
reaction of spreading and deformations of states on expectation values which are essential 
for dynamical effective equations [12]. An additional aspect arises for generally covariant 
situations where not all variables can be peaked in a semiclassical state as it would be the 
case in a common kinematical coherent state. Some of the phase space variables will have 
to play the role of internal "clocks" in which evolution of expectation values as well as 
quantum variables is measured. Thus, when constructing kinematical coherent states to 
check the classical limit, they must not be peaked on all phase space variables; a choice of 
clock has to be made before the calculation. Then, quantum variables also back-react on 
the change of the clock. 

Often, investigations of classical limits based on kinematical coherent states are moti- 
vated by well-known constructions of the harmonic oscillator or free quantum field theories. 
The behavior of quantum variables or of dynamical coherent states in general can, how- 
ever, be very different from the well-studied aspects of the free systems. Such systems or 
small deviations from them with anharmonicity can well be studied by coherent state tech- 
niques. But gravity is very different and not expanded around a set of harmonic oscillators. 
In fact, gravity with its unbounded Hamiltonian even lacks a ground state or vacuum to 
expand around. The bounce model solved in [20], for instance, shows that the spreads 
change exponentially rather than being constant or at least periodic as it happens for the 
harmonic oscillator. The resulting semiclassical picture is very different from that provided 
by harmonic oscillator coherent states. This must be taken into account in semiclassical 
analyses; effective theory provides suitable means to study such situations in sufficiently 
general terms as initiated in this paper. 
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5.3 Collective graviton 

The constructions indicate a picture of the classical limit of quantum gravity where linear 
metric modes appearing in the evolution equations are not basic excitations of a quantum 
field. They arise, rather, as collective excitations out of the underlying discrete quantum 
theory. At a basic level, degrees of freedom are encoded in quantum numbers fiyj while 
the scalar mode, for instance, is obtained through the difference between such a local label 
and the average value on the whole lattice. The classical modes thus arise as non-local, 
effective excitations out of the underlying quantum state ^. This shows in a well-defined 
sense how classical degrees of freedom are obtained as collective excitations, analogously to 
phonons in a crystal. That the correct classical dynamics results for these collective modes 
is demonstrated, for instance, by the derivation of Newton's potential for perturbations on 
a flat isotropic background in [6]. 

6 Summary 

Together with [121 [Zl E] we have shown in this paper that techniques are now available 
to derive effective equations of cosmological perturbation theory. The geometrical back- 
ground on which cosmological perturbations are defined is introduced through a class of 
states, rather than being used to set up the quantization. Background independent quan- 
tum properties thus remain, but one can make use of perturbation expansions for explicit 
calculations. The role of quantum labels and basic variables is clear from this procedure, 
which determines the type as well as order of magnitudes of correction terms which re- 
mained obscure previously. The inhomogeneous treatment including all relevant modes 
allows us to see all possible correction terms. Note, for instance, that since the isotropic 
expression for the co-triad is finite on small scales even classically, it would not contribute 
a correction function to the gravitational part of the Hamiltonian constraint in a purely 
isotropic setting. When isotropic expressions are quantized directly, a cancellation of the 
inverse isotropic triad component hides possible quantum effects of the full constraint. This 
does not arise in our context starting from an inhomogeneous lattice quantization. Thus, 
complete corrections are obtained in reliable form. 

Many different regimes are still to be explored to obtain a full overview of all effects. 
Moreover, gauge issues have to be investigated which is of relevance for the full quantum 
theory, too. While general effective equations including the relevant inhomogeneities are 
now available and orders of correction terms can be estimated, one still has to use them with 
care since they are not yet formulated for gauge invariant perturbations. In this context 
one should notice that not only evolution equations but also gauge transformations are 
determined by the constraints and thus modified by quantum corrections. It is thus not 
possible to use classical expressions for gauge invariant quantities since they will receive 
additional corrections. These issues are currently being studied to complete the derivation 
of equations with quantum corrections. The strategy for computing those corrections from 
quantum operators has been provided in this paper. 
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When speaking of quantum corrections to classical equations it is clear that to zeroth 
order the correct classical limit has to be satisfied. In fact, what we have shown here implies 
that loop quantum gravity has the correct classical limit for scalar modes in longitudinal 
gauge propagating on a spatially fiat background. In the process, we have demonstrated 
which steps must be involved in such a detailed calculation, most importantly a continuum 
limit but also a slowly varying field approximation. Extensions to other modes and gauges, 
and different backgrounds, can be done by the same techniques but are technically more 
involved to do explicitly. Nevertheless, it is clear from the construction that the correct 
classical limits will also be reproduced in those cases. More precisely, we have shown 
in Sec. 15.2.21 that there are always ranges of the basic lattice variables such that quantum 
corrections are small in nearly classical situations of low energy density and small curvature. 
In more energetic cosmological situations, those ranges can shrink to narrow intervals such 
that significant quantum corrections can be expected [6]. 

This demonstration of the correct classical limit crucially rests on a new understanding 
of effective theory [12] . Although it has not yet been formulated fully for field theories (but 
see [H]), this scheme is applicable here due to the ultraviolet cut-off of quantum gravity. 
On any lattice state we have only finitely many degrees of freedom in any compact spatial 
volume to which the quantum mechanical techniques of [12] directly apply. While loop 
quantum gravity does not possess a sharp cut-off but is rather based on arbitrary graphs 
in space any of which can occur in a general superposition [37], effective equations are 
always defined with respect to a single class of states. The physical state thus determines 
the cut-off dynamically. We have certainly not used explicit physical solutions of the 
Hamiltonian constraint but rather computed effective equations from general states. If a 
physical solution were available, the labels would be determined explicitly and fix the 
order of correction terms completely. Moreover, a full solution would determine how the 
lattice itself changes by the creation of new vertices in terms of an internal clock such as the 
total volume. This graph- changing nature seems to be one of the most important effects 
to be understood especially for late-time evolution in cosmology, or any dynamical issue 
relevant for large spatial slices. Although such a full solution seems currently out of reach, 
models and effective analyses already provide quite detailed information on the dynamics 
of background independent quantum gravity in cosmologically relevant regimes. 
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